{smcl}
{com}{sf}{ul off}{txt}{.-}
       log:  {res}E:\partiesdata\jop02rep.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res} 2 Aug 2002, 16:01:06
{txt}
{com}. *This log file replicates the results of Dan Reiter and Erik R. Tillman, "Public, Legislative, and Executive
. *Constraints on the Democratic Initiation of Conflict," Journal of Politics 64 (August 2002): 810-826.
. *Analysis was conducted in STATA 7.0.  Questions can be directed to Dan Reiter at dreiter@emory.edu.
. *This log file will replicate the results of Table 1 on p. 821.  Full description of research design and
. *variables is of course available in the article.
. use "E:\partiesdata\jop02rep.dta", clear
{txt}
{com}. *First, Model 1 of Table 1.
. xtgee    midindm majpow autbor polity if year>1918 & year<1993, fam(bin) link(probit) corr(ar1) t(year) force robust

{txt}Iteration 1: tolerance = {res}.02385192
{txt}Iteration 2: tolerance = {res}.00175295
{txt}Iteration 3: tolerance = {res}.000099
{txt}Iteration 4: tolerance = {res}6.001e-06
{txt}Iteration 5: tolerance = {res}3.692e-07

{txt}GEE population-averaged model                   Number of obs      = {res}     1721
{txt}Group and time vars:            {res}ccode year      {txt}Number of groups   = {res}       35
{txt}Link:                               {res}probit      {txt}Obs per group: min = {res}        3
{txt}Family:                           {res}binomial                     {txt}avg = {res}     49.2
{txt}Correlation:                         {res}AR(1)                     {txt}max = {res}       74
                                                {txt}Wald chi2({res}3{txt})       = {res}    55.72
{txt}Scale parameter:                 {res}        1      {txt}Prob > chi2        =    {res}0.0000

                            {txt}(standard errors adjusted for clustering on ccode)
{hline 13}{c TT}{hline 64}
             {c |}             Semi-robust
     midindm {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
      majpow {c |}  {res} 1.446246   .2006251     7.21   0.000     1.053028    1.839464
      {txt}autbor {c |}  {res} .2523288    .101682     2.48   0.013     .0530359    .4516218
      {txt}polity {c |}  {res} .0580762    .058216     1.00   0.318    -.0560251    .1721774
       {txt}_cons {c |}  {res}-2.026109    .587303    -3.45   0.001    -3.177201   -.8750159
{txt}{hline 13}{c BT}{hline 64}

{com}. *Now, Model 2.
. xtgee  midindm  vanpart   eff_ptie  majdum coaldum presdum prempres assind treatrat majpow autbor if year>1918 & year<1993, fam(bin) link(probit) corr(ar1) t(year) force robust

{txt}Iteration 1: tolerance = {res}.1398609
{txt}Iteration 2: tolerance = {res}.02537971
{txt}Iteration 3: tolerance = {res}.00109677
{txt}Iteration 4: tolerance = {res}.00014404
{txt}Iteration 5: tolerance = {res}.00001524
{txt}Iteration 6: tolerance = {res}7.371e-07

{txt}GEE population-averaged model                   Number of obs      = {res}     1363
{txt}Group and time vars:            {res}ccode year      {txt}Number of groups   = {res}       28
{txt}Link:                               {res}probit      {txt}Obs per group: min = {res}        3
{txt}Family:                           {res}binomial                     {txt}avg = {res}     48.7
{txt}Correlation:                         {res}AR(1)                     {txt}max = {res}       74
                                                {txt}Wald chi2({res}9{txt})       = {res}   207.49
{txt}Scale parameter:                 {res}        1      {txt}Prob > chi2        =    {res}0.0000

                            {txt}(standard errors adjusted for clustering on ccode)
{hline 13}{c TT}{hline 64}
             {c |}             Semi-robust
     midindm {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
     vanpart {c |}  {res}-.0126103   .0053703    -2.35   0.019    -.0231358   -.0020847
    {txt}eff_ptie {c |}  {res} -.002886   .0906549    -0.03   0.975    -.1805664    .1747944
      {txt}majdum {c |}  {res} .1736322   .2236612     0.78   0.438    -.2647357    .6120001
     {txt}coaldum {c |}  {res} .2017317   .2397001     0.84   0.400     -.268072    .6715353
     {txt}presdum {c |}  {res} .0220547   .4062441     0.05   0.957    -.7741692    .8182786
    {txt}prempres {c |}  {res}-.4880111   .2245778    -2.17   0.030    -.9281756   -.0478467
      {txt}assind {c |}  {res}-1.126466   .2697834    -4.18   0.000    -1.655232   -.5977004
    {txt}treatrat {c |}  {res}-.3045961   .1639302    -1.86   0.063    -.6258934    .0167011
      {txt}majpow {c |}  {res}  1.54259   .2583133     5.97   0.000     1.036306    2.048875
      {txt}autbor {c |}  {res} .2826882   .0778863     3.63   0.000     .1300339    .4353425
       {txt}_cons {c |}  {res}-.7026606   .4842636    -1.45   0.147      -1.6518    .2464787
{txt}{hline 13}{c BT}{hline 64}

{com}. *Next, Model 3.
. xtgee  midindm  vanpart   eff_ptie  majdum coaldum presdum  semipres treatrat majpow autbor if year>1918 & year<1993, fam(bin) link(probit) corr(ar1) t(year) force robust

{txt}Iteration 1: tolerance = {res}.14335111
{txt}Iteration 2: tolerance = {res}.02024109
{txt}Iteration 3: tolerance = {res}.00072588
{txt}Iteration 4: tolerance = {res}.00011425
{txt}Iteration 5: tolerance = {res}.00001115
{txt}Iteration 6: tolerance = {res}7.345e-07

{txt}GEE population-averaged model                   Number of obs      = {res}     1363
{txt}Group and time vars:            {res}ccode year      {txt}Number of groups   = {res}       28
{txt}Link:                               {res}probit      {txt}Obs per group: min = {res}        3
{txt}Family:                           {res}binomial                     {txt}avg = {res}     48.7
{txt}Correlation:                         {res}AR(1)                     {txt}max = {res}       74
                                                {txt}Wald chi2({res}9{txt})       = {res}   111.58
{txt}Scale parameter:                 {res}        1      {txt}Prob > chi2        =    {res}0.0000

                            {txt}(standard errors adjusted for clustering on ccode)
{hline 13}{c TT}{hline 64}
             {c |}             Semi-robust
     midindm {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
     vanpart {c |}  {res}-.0104474   .0050436    -2.07   0.038    -.0203326   -.0005621
    {txt}eff_ptie {c |}  {res}-.0151781   .0904252    -0.17   0.867    -.1924081     .162052
      {txt}majdum {c |}  {res} .1493097   .2221978     0.67   0.502      -.28619    .5848095
     {txt}coaldum {c |}  {res} .1854494   .2406672     0.77   0.441    -.2862497    .6571485
     {txt}presdum {c |}  {res} .0279669   .4070491     0.07   0.945    -.7698348    .8257685
    {txt}semipres {c |}  {res}-.6020274   .2360841    -2.55   0.011    -1.064744   -.1393111
    {txt}treatrat {c |}  {res}-.3032556   .1641821    -1.85   0.065    -.6250467    .0185354
      {txt}majpow {c |}  {res} 1.572963    .265868     5.92   0.000     1.051871    2.094055
      {txt}autbor {c |}  {res} .2874745   .0773398     3.72   0.000     .1358913    .4390578
       {txt}_cons {c |}  {res}-.7605716   .5047928    -1.51   0.132    -1.749947    .2288041
{txt}{hline 13}{c BT}{hline 64}

{com}. *Lastly, Model 4.

{com}. xtgee   mnin45dm vanpart   eff_ptie  majdum coaldum presdum prempres assind treatrat majpow autbor if year>1918 & year<1993, fam(bin) link(probit) corr(ar1) t(year) force robust

{txt}Iteration 1: tolerance = {res}.21238074
{txt}Iteration 2: tolerance = {res}.06032816
{txt}Iteration 3: tolerance = {res}.00832424
{txt}Iteration 4: tolerance = {res}.00155259
{txt}Iteration 5: tolerance = {res}.00020903
{txt}Iteration 6: tolerance = {res}.00003918
{txt}Iteration 7: tolerance = {res}3.703e-06
{txt}Iteration 8: tolerance = {res}6.062e-07

{txt}GEE population-averaged model                   Number of obs      = {res}     1363
{txt}Group and time vars:            {res}ccode year      {txt}Number of groups   = {res}       28
{txt}Link:                               {res}probit      {txt}Obs per group: min = {res}        3
{txt}Family:                           {res}binomial                     {txt}avg = {res}     48.7
{txt}Correlation:                         {res}AR(1)                     {txt}max = {res}       74
                                                {txt}Wald chi2({res}9{txt})       = {res}   143.22
{txt}Scale parameter:                 {res}        1      {txt}Prob > chi2        =    {res}0.0000

                            {txt}(standard errors adjusted for clustering on ccode)
{hline 13}{c TT}{hline 64}
             {c |}             Semi-robust
    mnin45dm {c |}      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
{hline 13}{c +}{hline 64}
     vanpart {c |}  {res}-.0120187   .0055327    -2.17   0.030    -.0228626   -.0011748
    {txt}eff_ptie {c |}  {res}-.0014211    .103268    -0.01   0.989    -.2038226    .2009803
      {txt}majdum {c |}  {res} .0449379   .1908869     0.24   0.814    -.3291935    .4190693
     {txt}coaldum {c |}  {res} .1843955   .1858587     0.99   0.321    -.1798808    .5486719
     {txt}presdum {c |}  {res} .0589727   .4312658     0.14   0.891    -.7862927    .9042381
    {txt}prempres {c |}  {res}-.3559427    .255747    -1.39   0.164    -.8571975    .1453122
      {txt}assind {c |}  {res}-1.168592   .3137564    -3.72   0.000    -1.783543   -.5536409
    {txt}treatrat {c |}  {res} -.302518   .2011354    -1.50   0.133    -.6967362    .0917002
      {txt}majpow {c |}  {res} 1.182541    .240246     4.92   0.000     .7116673    1.653414
      {txt}autbor {c |}  {res} .2088073   .0901546     2.32   0.021     .0321075     .385507
       {txt}_cons {c |}  {res}-.7934124   .5264754    -1.51   0.132    -1.825285    .2384604
{txt}{hline 13}{c BT}{hline 64}

{com}. log close
       {txt}log:  {res}E:\partiesdata\jop02rep.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res} 2 Aug 2002, 16:07:15
{txt}{.-}
{smcl}
{txt}{sf}{ul off}
